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We propose using the optomechanical interaction to create artificial magnetic fields for photons on 
a lattice. The ingredients required are an optomechanical crystal, i.e. a piece of dielectric with the 
right pattern of holes, and two laser beams with the right pattern of phases. One of the two proposed 
schemes is based on optomechanical modulation of the links between optical modes, while the other 
is a lattice extension of optomechanical wavelength-conversion setups. We illustrate the resulting 
optical spectrum, photon transport in the presence of an artificial Lorentz force, edge states, and 
the photonic Aharonov-Bohm effect. Moreover, we briefly describe the gauge fields acting on the 
synthetic dimension related to the phonon/photon degree of freedom. 


Light interacting with nano-mechanical motion via the 
radiation pressure force is studied in the field of optome- 
chanics. The field has seen rapid progress in the last few 
years (see the recent review pQ). So far, most experimen¬ 
tal achievements have been realized in setups comprising 
one optical mode coupled to one vibrational mode. Ob¬ 
viously, one of the next frontiers will be the combination 
of many such optomechanical cells into an optomechan¬ 
ical array, enabling the optical in-situ investigation of 
(quantum) many-body dynamics of interacting photons 
and phonons. Many experimental platforms allow to be 
scaled up to arrays. However, optomechanical crystals 
seem to be the best suited candidate at the present stage. 
Optomechanical crystals are formed by the periodic spa¬ 
tial patterning of regular dielectric and elastic materi¬ 
als, resulting in an enhanced coupling between optical 
and acoustic waves via moving boundary or electrostric- 
tion radiation pressure effects. Two-dimensional (2D) op¬ 
tomechanical crystals with both photonic and phononic 
bandgaps [2j can be fabricated by standard microfabri¬ 
cation techniques through the lithographic patterning, 
plasma etching, and release of a thin-film material [3]. 
These 2D crystals for light and sound can be used to cre¬ 
ate a circuit architecture for the routing and localization 
of photons and phonons [3HZ]. 

Optomechanical arrays promise to be a versatile plat¬ 
form for exploring optomechanical many-body physics. 
Several aspects have already been investigated theoret¬ 
ically, e.g. synchronization [SHIS, long-range interac¬ 
tions m na , reservoir engineering m, entanglement 
mi na, correlated quantum many-body states m, slow 
light m, transport in a ID chain mi, and graphene-like 
Dirac physics m- ■ 

One of the central aims in photonics is to build waveg¬ 
uides that are robust against disorder and do not display 
backscattering. Recently there have been several pro¬ 
posals fl9lf23j to engineer non-reciprocal transport for 
photons. On the lattice, this corresponds to an arti¬ 
ficial magnetic field,which would (among other effects) 
enable chiral edge states that display the desired robust¬ 
ness against disorder. First experiments have shown such 
edge states [24h26] . These developments in photonics are 


related to a growing effort across various fields to produce 
synthetic gauge fields for neutral particles [271429] , 

In this paper we will propose two schemes to gener¬ 
ate an artificial magnetic field for photons on a lattice. 
In contrast to any previous proposals or experiments for 
photonic magnetic fields on a lattice, these would be con¬ 
trolled all-optically and, crucially, they would be tunable 
in-situ by changing the properties of a laser field (fre¬ 
quency, intensity, and phase pattern). They require no 
more than a patterned dielectric slab illuminated by two 
laser beams with suitably engineered optical phase fields. 
The crucial ingredient is the optomechanical interaction. 

On the classical level, a charged particle subject to a 
magnetic field experiences a Lorentz force. In the quan¬ 
tum regime, the appearance of Landau levels leads to the 
integer and fractional quantum Hall effects, where topo¬ 
logically protected chiral edge states are responsible for a 
quantized Hall conductance. On a closed orbit, a particle 
with charge q will pick up a phase that is given by the 
magnetic flux <f> through the circumscribed area, where 
= (q/h) fB-dS in units of the flux quantum, with 
B denoting the magnetic field. On a lattice, a charged 
particle hopping from site i to j acquires a Peierls phase 
(fij = (q/K) Adr determined by the vector potential 
A. Conversely, if we can engineer a Hamiltonian for neu¬ 
tral particles containing arbitrary Peierls phases, 

Hhop = hJ ^ e^otoi + h.c., (1) 

(U) 

we are able to produce a synthetic magnetic field. Here 
di is the (bosonic) annihilation operator on lattice site 
i. We note in passing that different phase configurations 
can lead to identical flux patterns, reflecting the gauge 
invariance of Maxwell’s equations under the transforma¬ 
tion A —>■ A + V£(r) for any scalar function £. 

Every defect in an optomechanical crystal mz] sup¬ 
ports a localized vibrational (annihilation operator b , 
eigenfrequency Do) and optical mode (a, frequency u; cav ) 
that interact via radiation pressure, giving rise to the 
standard optomechanical interaction [T]: 

Hint = -hg 0 a ] a(b ] + b). 
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Figure 1: Proposed schemes to create a photonic gauge field in 
optomechanical arrays by engineering photon hopping phases, 
(a) Modulated link scheme, (b) Corresponding optical spec¬ 
trum of a row with relevant sidebands (dashed). Driven vi¬ 
brational modes (yellow) optomechanically modulate the fre¬ 
quency of optical link modes. Tunneling photons are thus 
up-converted to the first sideband and pick up the phase of 
the modulation. Arrows in (b) indicate the resonant pho¬ 
ton transmission process in a row. (c) Wavelength conversion 
scheme and (d) corresponding optical spectrum: Neighboring 
modes in a row couple to a vibrational mode (yellow) optome¬ 
chanically (red lines, denoting the linearized optomechanical 
interaction). Two lasers, driving the optical modes close to 
the red sidebands (wiggly arrows in d), give rise to resonant 
photon-phonon-photon conversion with the phase set by the 
lasers. Different rows are connected via simple photon hop¬ 
ping (blue lines), without a phase, in both schemes (a,c). The 
indicated phase configuration corresponds to a constant mag¬ 
netic field. 


This can be utilized in two basic ways to introduce phases 
for the hopping of photons. First, one can drive the 
optical mode by a control laser (frequency ujl) close to 
the red sideband, ujl ~ — ^o- Following the stan¬ 

dard procedure of linearization and rotating wave ap¬ 
proximation (RWA) [I] one recovers a swap Hamiltonian, 
gam + h.c., in which the phase (j) of the coupling g E C 
is set by the control laser phase. We will show below 
how this can be used to create a photonic gauge field. 
There is, however, also a second route, namely, driv¬ 
ing the vibrational mode into a large amplitude coherent 
state, (b(t)) = |/3|e“*( fit+ ^), using the radiation pressure 
force. These oscillations then weakly modulate the opti¬ 
cal eigenfrequency, cj cav (t) = u; caV;0 + 2g 0 \(}\ cos (Sit + 0), 
with the phase </> set by the oscillations. Again, in a suit¬ 
able setting this will lead to an artificial magnetic field 
for the photons. We now describe both methods in turn. 

Modulated link scheme. - Recently Fang et. al. 
[la mi proposed to create a photonic gauge field by 
electro-optically modulating the photon hopping rate 
Jij = J cos (Sit + between neighboring cavities. This 
would require locally wired electrodes for each link of 


the lattice. Here we propose a potentially more power¬ 
ful all-optical implementation of that idea. We employ 
optomechanically driven photon transitions, as first dis¬ 
cussed in m, but extended to a scheme with modulated 
interface modes, depicted in Figure]]] (a). We now discuss 
the leftmost three optical modes in the first row, a a, &/, 
clb (from left to right), exemplary of the full grid. Their 
coherent dynamics is governed by the Hamiltonian 

H/h = Uid\di + uji{t)a\di — J(a\aA + h.c.). 

i=A,B 

(3) 

The terms describe, in this order, the first (A) and third 
(B) optical mode, the temporally modulated interface 
mode (I), and its tight binding coupling to the neigh¬ 
boring A and B modes with photon tunneling rate J. 
As discussed further below, the eigenfrequency uj of the 
interface mode should be well separated from the eigen- 
frequencies of the adjacent A and B modes, for the transi¬ 
tion A-l-B to be virtual. The interface mode is optome¬ 
chanically coupled to a mechanical mode, which itself is 
driven into a large amplitude coherent state. As men¬ 
tioned above, this gives rise to a weak modulation of its 
optical eigenfrequency, uji(t) = a)/ + 2go\/3\ cos (Sit + </>), 
with the phase (j) set by the driving. The required me¬ 
chanical driving is easily generated by two-tone laser ex¬ 
citation at a frequency difference SI. The beating be¬ 
tween the laser beams gives rise to a sinusoidal radia¬ 
tion pressure force, which drives the mechanical mode. 
If ujb = + SI, then a photon hopping from site A 

to B picks up the phase </> of the modulation: Start¬ 
ing from a^, it tunnels into dj where it is inelasti- 
cally up-scattered into the first sideband by the modu¬ 
lation and subsequently tunnels into a# resonantly, as 
shown by the spectrum in Figure [ljb). We can de¬ 
rive an effective Hamiltonian, hJ^e^a^aA + h.c. for 
this process by integrating out the interface mode I us¬ 
ing Floquet perturbative methods to third order (see 
Appendix \M. For the effective hopping rate we find 
Jeff = go\/3 \J 2 /[(ua - w;)(wb - w/)], to leading order 
in J and \/3\. Concatenating such three-mode blocks, we 
create a linear chain (the first row in Figure |TJl) , with its 
optical spectrum schematically depicted in Figure [ljb). 
Every time a photon hops to the right, it is up-converted 
and picks up the phase of the drive. To obtain a 2d grid, 
we stack identical chains and connect neighboring rows 
by direct photon hopping (whose rate must be chosen 
to equal J e ff, to obtain isotropic hopping), as depicted 
in Figure ]ljb). The phase configuration in Figure B a ) 
corresponds to a constant magnetic field. Note that in 
contrast to the general Hamiltonian ([!]), this scheme does 
not allow for phases when hopping between rows, yet it 
is still possible to achieve an arbitrary flux through every 
plaquette. Hence, arbitrary magnetic fields can be gen¬ 
erated, provided one can control the driving laser phase 
at every interface mode. With the help of wave front en¬ 
gineering, this can be achieved with no more than two 
lasers: A homogeneous ’carrier’ beam E% = Eioe~ luJLt 
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Figure 2 : Comparison between the experimentally accessible optical spectrum (LDOS) of the ideal effective Hofstadter model 
(a) and the actual results from the proposed modulated link scheme (b-d), for different optical decay rates and magnetic fluxes. 
The simulation results indicate that the scheme works even beyond the perturbative regime. The resulting Hofstadter butterfly 
could be observed by a local tapered fiber probe. Modulation of links produces higher side bands (d). The phase configuration 
corresponds to a constant magnetic flux <f> per plaquette, see Fig. [l] [Parameters: grid: 10 x 10 (a), 12 x 12 (b-d); J e ff = O.IO 8 D 0 
(a-d); ft — O.OlDo (a,c,d), ft = 0.05Do (b); J — 0.3Do (b-d); go|/3| = 0.3Do (b-d); optical eigenfrequencies (relative to first mode) 
in a row (left to right) including interface modes: (0, 0.5,1,1.5,... )Qo (b-d).] 


and a ’modulation’ beam E 2 = E 2 oe~^ CJL+n ^ t ~ l ^ x,y \ 
with an imprinted phase pattern (j)(x,y). Interfer¬ 
ence yields the desired temporally modulated intensity 
\E 10 f + |£ 2 o | 2 + 2Re[Et 0 E 20 e- i ^ t +^ x ^], exerting a 
radiation force with a site-dependent phase. Care has 
to be taken to avoid exciting other vibrational modes 
(those not at the interface mode), by engineering them 
to have different mechanical frequencies. To this end, 
the driving frequency D would usually be chosen close 
to the mechanical eigenfrequency Do, so the mechanical 
amplitude is enhanced by the mechanical quality factor 
and is thus much larger than any spurious amplitude in 
other (off-resonant) modes. By engineering the intensity 
pattern \E 2 o(x,y)\ as well, one could suppress any such 
unwanted effects even further. 


Wavelength conversion scheme. - There is another, 
alternative way of engineering an optical Peierls phase, 
and it is related to optomechanical wavelength conver¬ 
sion [32j|33]. In wavelength conversion setups, low fre¬ 
quency photons in one mode are up-converted to a higher 
frequency in another mode by exploiting the modes’ mu¬ 
tual optomechanical coupling to a vibrational mode. We 
propose to scale up this idea into a grid as depicted in 
Figure [lie). The leftmost three modes in the first row 
depict (in this order) an optical mode (annihilation op¬ 
erator clai frequency lja), a mechanical mode ( 6 , Do) 
and another optical mode (a#, ujb 7 ^ ^a)- The me¬ 
chanical mode couples optomechanically to both opti¬ 
cal modes. A and B are driven by a laser with fre¬ 
quency ujli and ul2, respectively: For mode A , we re¬ 
quire cja ~ u li = Do + £ = D where ujli denotes the 
driving laser’s frequency and S <C fto is the detuning 
from the red sideband. For mode H, a similar relation 
mb — ul 2 = ft holds, as depicted in the spectrum in Fig¬ 
ure Bd). After application of the standard linearization 
and RWA procedure [l], the dynamics in a frame rotating 


with the drive is governed by the Hamiltonian 

H/h = Cl ajcii + QoPb — (g\aAb^ + gB^ B b + h.c.). 

i=A,B 

(4) 

Elimination of the mechanical mode leads to an effec¬ 
tive Hamiltonian hJ^e 1 ^c^ b cla + h.c. to leading order in 
\9a,b\/S, with effective hopping rate J eff = \qa\\9b\/S 
and hopping phase 4> = 4>b ~ 4>a • Here, <pA and <j)B are 
the phases of the linearized optomechanical interaction, 
of the form gA = \gA\E^ A : which are set by the phase of 
the laser drive at the corresponding site. Connecting al¬ 
ternating A and B sites by mechanical link modes yields 
a row whose spectrum is depicted in Figure [T] (d). As 
in the previous scheme, we can simply connect rows by 
photonic hopping without phases (at a rate J e ff) to yield 
a 2D grid. Phase front engineering of the two driving 
lasers is sufficient to realize arbitrary magnetic fields for 
photons in the grid. We note that the scheme also works 
for driving far away from the red sideband (yielding en¬ 
hanced values of D and thereby J e fr; see below), though 
that requires stronger driving. 

Another optomechanical scheme for non-reciprocal 
photon transport that could potentially be extended to 
a lattice is based on optical microring resonators |23j, 
but the connection of these rather large rings via waveg¬ 
uides would presumably result in a more complicated and 
less compact structure than what can be done with the 
photonic-crystal based approaches analyzed here. 

We now discuss the limitations imposed on the achiev¬ 
able effective hopping J e ^. The important end result will 
be that J e ff is limited to about the mechanical frequency 
Dq, even though perturbation theory would seem to im¬ 
ply a far smaller limit (for possible technical limitations 
connected to the driving strength, see the Supplementary 
Information). 

We denote as e <C 1 the order of the three small 
parameters J/\uja ~ ^ 1 1, J/\^b — cj/|, and go\/3\/Q in 
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the modulated link scheme (Fig. [l] a,b). Then the 
effective coupling strength in the perturbative regime 
reads J e q = 0{e 3 )Vt. Even though the modulation fre¬ 
quency D need not equal the eigenfrequency Do, they 
should usually be close to yield a significant mechani¬ 
cal response and avoid other resonances. For the wave¬ 
length conversion scheme, where \gA,B\/$ = 0(e), we 
recover J e q = 0(e 2 )S = 0(e 3 )Do, since RWA requires 
5/f2o to be small as well. In any experimental realiza¬ 
tion, photons will decay at the rate a. Thus they travel 
rsj Jeff / ^ rsj (Do/a) 0(e 3 ) sites. In order for the photons 
to feel the magnetic field (or to find nontrivial transport 
at all), this number should be larger than 1. That pre¬ 
cludes being in the deep perturbative limit e 1, even 
for a fairly well sideband-resolved system (where typi¬ 
cally a ~ 0.1 Do). Similar considerations apply for other 
proposed (non-optomechanical) schemes based on mod¬ 
ulation [22] . 

We now explore numerically the full dynamics, beyond 
the perturbative limit. The optical local density of states 
(LDOS) is experimentally accessible by measuring the 
reflection when probing an optical defect mode via a ta¬ 
pered fiber, and it reveals the spectrum of the Hamilto¬ 
nian. It thus provides a reasonable way to asses the valid¬ 
ity of the effective Hamiltonian beyond the perturbative 
limit. Figure [2] (a) shows the LDOS in the bulk calcu¬ 
lated with the ideal effective Hofstadter model 0 for a 
spatially constant magnetic field, depicting the famous 
fractal Hofstadter butterfly structure m- For compar¬ 
ison, we plot the LDOS of the modulated link scheme 
in Figure [2|b,c). It is obtained by calculating numeri¬ 
cally the Floquet Green’s function of the full equations 
of motion (with time-periodic coefficients), see Appendix 
iBl The results indicate that the scheme works even for 
Jeff ~ 0.1D > a, although perturbation theory clearly 
breaks down in this regime. We stress that the butterfly 
in Fig. [2|b, c) could even be observed experimentally at 
room temperature, since the spectrum is insensitive to 
thermal fluctuations. One would also observe sidebands, 
see Figure [2jd) . Similar results hold for the wavelength- 
conversion scheme (not shown here). 

In addition to measuring the optical spectrum, it is 
also possible to look at photon transport in a spatially 
resolved manner, by injecting a probe laser locally and 
then imaging the photons leaving the sample. This pro¬ 
vides another way to observe the effects of the artificial 
gauge field, which gives rise to distinct transport phe¬ 
nomena as depicted in Figure [3^ a,b). For small magnetic 
fields, |0| <C 27r, the dynamics can be understood in the 
continuum limit when probing the bulk: One recovers the 
standard Landau level picture for electrons in a constant 
magnetic field [34, 135], with effective mass m* = h/2Ja 2 
and cyclotron frequency u; cyc = 20 J, where a is the lat¬ 
tice constant. In Figure [3] (a) the n = 1 Landau level is 
selected via the probe’s detuning A p with respect to the 
drive. The circles indicate the semi-cla ssical cyclotron 
orbits with radius R = ay (2 n + 1)/0. In this semi- 
classical picture, the momentum of a photon injected 



Figure 3: Microscopic simulation of the wavelength conversion 
scheme, Eq. ([4]), indicating its feasibility: Spatial distribution 
of light intensity upon local injection of a probe laser in the 
bulk (a) and at the edge (b), for a constant artificial magnetic 
field. Bulk transport (a) is governed by Landau levels and can 
be understood as a superposition of classical cyclotron orbits 
(yellow circles) for different momentum directions, (b) At 
the edges robust edge channels exist, (c) Optical Aharonov- 
Bohm effect in minimal symmetric setup: (d) The interfer¬ 
ence pattern (normalized probe laser transmission intensity) 
is shifted by the magnetic flux through the ring. [Parameters: 
22x22 grid (a,b); S = 0.3D o (a,b), S = 0.1D o (d); g = 0.2D o 
(a,b), g = O.OlHo (d); a = 0.01D o ; r = a/10; 4> = 2tt/8 
(a,b); J = 0.13D o (a,b), J = 0.001D o (d), A p = 1.278D 0 (a), 
A p = 1.260D o (b), Ap = 1.103D o (d)] 


locally at a site in the bulk is equally distributed over 
all directions, since the position is well-defined. Thus, 
the observed response resembles a superposition of semi- 
classical circular Lorentz trajectories with different initial 
velocity directions. A probe injected closer to the edge 
excites chiral integer Quantum Hall Effect edge states, 
see Fig. [3^b). 

The Aharonov-Bohm effect [36] is one of the most in¬ 
triguing features of quantum mechanics. In an interfer¬ 
ometer, electrons can acquire a phase difference deter¬ 
mined by the magnetic flux enclosed by the interfering 
pathways, even though they never feel any force due to 
the magnetic field. Figure |3][c) depicts a setup that is 
based on the wavelength conversion scheme and realizes 
an optical analog of the Aharonov-Bohm effect: A lo¬ 
cal probe is transmitted via two pathways, leading to an 
interference pattern in the transmission. The pattern is 
shifted according to the flux through the ’ring’, see Fig. [3] 
(d), confirming the effect. 

All the effects displayed in Fig. [3] have been simulated 
numerically for the wavelength conversion scheme, see 
Appendix [Cj but similar results hold for the modulated- 
link scheme. 

So far we have analyzed schemes to engineer hopping 
phases for photons. We now ask about situations where 
the phonons are not only employed as auxiliary virtual 
excitations, but rather occur as real excitations, which 
can be interconverted with the photons. This means, in 
addition to the modes making up the lattices described 
above (in either of the two schemes), we now consider 
on-site vibrational modes bj coupled optomechanically to 
the corresponding optical modes aj. Using the standard 
approach p], we arrive at a linearized optomechanical 
interaction of the form —ga^bj + h.c.. Moreover, to be 
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Figure 4: Optomechanical gauge fields within the concept 
of synthetic dimensions, (a) The optomechanical coupling, g , 
can be viewed as connecting sites along a synthetic dimension 
(photons vs. phonons). A phase for the photon hopping, 
engineered using the schemes from above, creates a flux in 
the optical plaquette (blue, top) and in the adjacent synthetic 
plaquette (gray). Hence, the magnetic field (black lines) in the 
full space is divergence-free, (b) Engineering exclusively the 
phases of g allows to create magnetic fields/fluxes, but only 
perpendicular to the synthetic dimension. Shining a single 
tilted laser on a ID chain yields a synthetic optomechanical 
ladder system with constant synthetic magnetic flux, (c) 2D 
array, with the held inside the (physical) plane generated by 
an arbitrary laser phase pattern. 


general (and generate nontrivial features connected to the 
gauge held structure), we will assume the neighboring 
phonon modes may also be coupled, as described by a 
tight-binding Hamiltonian of the form —iF bjbi + 

h.c.. 

When discussing the effects of gauge helds in such a 
setting, the system is best understood within the con¬ 
cept of ’synthetic’ dimensions [33 EH]. The optomechan¬ 
ical interaction can be viewed in terms of an extension 
of the ID or 2D lattice into such an additional synthetic 
dimension. In our case, this dimension only has two dis¬ 
crete locations, corresponding to photons vs. phonons. 
In that picture, the optomechanical interaction, convert¬ 
ing photons to phonons, corresponds to a simple hopping 
between sites along the additional direction. Figure [4] (a) 
sketches this for an optomechanical ring: photons and 
phonons represent two layers separated along the syn¬ 
thetic dimension. Applying any of our two previously 
discussed schemes, a photon hopping from site i to j will 
acquire a phase <pij = dr A. The gauge field A must 
now be viewed as a vector field in this new 3d space, 
where one of the dimensions is synthetic. A finite hop¬ 
ping phase (p at one of the optical links creates a magnetic 
flux through the optical plaquette as desired, see Fig. [4] 
(a). However, and this is the important point, since the 
magnetic field B is divergence-free, the field must pen¬ 
etrate at least one additional plaquette, causing the op¬ 
posite magnetic flux in the synthetic dimension (assum¬ 
ing g G M). In general, realizing that there is this kind 
of behaviour is crucial to avoid puzzles about seeming 
violations of gauge symmetry in situations with photon 
magnetic fields in optomechanical arrays. It is necessary 


to keep track of the full vector potential in the space that 
includes the synthetic dimension. 

We now take a step back, getting rid of the previ¬ 
ously discussed engineered schemes that required two 
lasers and some arrangement of ’link’ modes. Rather 
we will consider simple optomechanical arrays, i.e. lat¬ 
tices of optical and vibrational modes, with photon and 
phonon tunnel coupling between modes and with the op¬ 
tomechanical interaction. We ask: What is the effect of 
an arbitrary, spatially varying optical phase field in the 
driving laser that sets the strength of the optomechanical 
coupling? It turns out that the resulting spatially vary¬ 
ing phase of the optomechanical coupling, gj = \gj\e l(pj , 
can be chosen to create arbitrary magnetic fields perpen¬ 
dicular to the synthetic dimension. A particularly sim¬ 
ple example is a simple linear chain of optomechanical 
cells. Shining a tilted laser (i.e. with a phase gradient, 
(fj = j -5(p) onto such a ID optomechanical array creates 
a constant magnetic flux through the plaquettes of the 
“optomechanical synthetic ladder” that can be drawn to 
understand the situation, cf. Figure [4] (b). The quantum 
mechanics of excitations tunneling between the two ’rails’ 
of the ladder (corresponding to photon-phonon conver¬ 
sion) is directly analogous to experiments on electron 
tunneling between parallel wires in a magnetic field m- 
The magnetic field shifts the momenta of the tunneling 
particles, giving rise to resonance phenomena when the 
shifted dispersion curves cj(k) of the excitations match. 
Via phase front engineering one could create arbitrary 
synthetic magnetic fields also in 2D grids, see Figure 
|4jc). We note, though, that this method is constrained 
since it cannot create directly magnetic fluxes through 
optical or mechanical plaquettes, and in general only the 
schemes discussed above provide full flexibility. On the 
other hand, if either the photon or phonon modes are 
occupied only virtually, then effective fluxes can still be 
generated for the remaining real excitations, even with a 
single laser, and this works best for phonons (see [40, 41]). 

We now discuss the most salient aspects of the ex¬ 
perimental realization. Both the ’butterfly’ optical spec¬ 
trum and spatially resolved transport can be probed us¬ 
ing homodyne techniques, which are insensitive to noise. 
Real-space optical imaging is feasible, as the defects are 
a few micrometers apart. The optical phase pattern can 
be engineered using spatial light modulators. No time- 
dependent changes of the pattern are needed here, since 
the time-dependence is generated via the beat-note be¬ 
tween the two laser beams. 

For the modulated-link scheme, the mechanical oscil¬ 
lation amplitude (3 used for the modulation should over¬ 
whelm any thermal fluctuations. In the example of Fig.[2j 
we assumed go(3 = 0.3Do. At recently achieved parame¬ 
ters [7] go/27r = 220 kHz and D 0 /27r = 9 GHz , this would 
imply (3 ~ 10 4 , i.e. a phonon number of 10 8 reached by 
driving, certainly larger than the thermal population. If 
we drive the mechanical vibration using a radiation pres¬ 
sure force oscillating at resonance (assuming the quoted 
go also for the optical mode used in that driving), then 
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we have /3 = 2gon c /Y , where n c is the circulating photon 
number and Y the mechanical damping rate. Given a 
mechanical quality factor of f^o/r = 2 • 10 5 , this requires 
n c ~ 10 3 photons for Fig. [ 2 J a realistic number. We note 
that thermal fluctuations of the mechanical amplitude 
give rise to a fractional deviation of y/fith /ft in T e ff, with 
a slow drift on the time scale T _1 . At typical temper¬ 
atures used in experiments, we have n t h 00 100 — 1000 , 
and so the fractional change is on the order of a percent, 
which will not noticeably impact transport. 

In the alternative wavelength-conversion scheme, one 
should strive for a large photon-enhanced optomechani¬ 
cal coupling rate g = g^a. A general estimate implies we 
always need the photon number to be larger than (k/ go) 2 
in order to see the butterfly spectrum and the transport 
effects. This condition (compatible with Fig. [3| would 
require a circulating photon number of around 3 • 10 5 for 
the parameters demonstrated in a recent successful wave¬ 
length conversion setup based on optomechanical crystals 
|32| . It is also important to estimate the unwanted in¬ 
flux of thermal excitations from the phonon subsystem 
into the photon subsystem, at least if the setup is to be 
applied in the quantum regime, for observing the trans¬ 
port of single photons in the presence of a magnetic field. 
In the wavelength-conversion scheme, there is a detun¬ 
ing 5 k between the red sideband of the laser and 
the phonon mode, such that photon-phonon conversion 
is suppressed. Nevertheless, it still happens at a rate yn, 
where 7 = g 2 /5 is the “cooling rate” (for the detuned case 
applicable here) and n is the number of phonons in the 
mode. Fortunately, this phonon number is also reduced 
by the very same off-resonant cooling process. Balanc¬ 
ing the inflow and outflow of excitations, we find that 
there will be a remaining unwanted photon occupation of 
^phot = ™thY /k due to the conversion of thermal phonons 
into photons, where n t h is the bulk thermal phonon occu¬ 
pation. The factor Y/k suppresses this number strongly, 
and it should be possible to reach the regime hp 1 J lot 1 
in low-temperature setups. 

Reducing fabrication-induced disorder will be crucial 
for any future applications of photonic crystals, includ¬ 
ing the one envisaged here (as well as other photonic 
magnetic field schemes). In first experimental attempts, 
the optical and mechanical disorder is on the percent 
level, which makes especially the fluctuations of the 
optical resonance frequencies significant. Nevertheless, 
strong reductions of the disorder will be possible by post¬ 
fabrication methods [421144] , such as local laser-induced 
oxidation. These are expected to reduce the fluctuations 
down to the level of 10“ 5 relative optical frequency fluc¬ 
tuations. This is enough to suppress the optical disorder 
to some fraction of the photon hopping rate J e ff ~ fl, 
which will enable near-ideal photon transport (e.g. An¬ 
derson localization lengths would be at least hundreds 
of sites, larger than the typical arrays). Disorder in the 
mechanical frequencies can be reduced by similar tech¬ 
niques, but is much less problematic, due to the difference 
in absolute frequency scales between optics and mechan¬ 


ics. 

Outlook - Optomechanical crystals represent an inter¬ 
esting system for the realization of artificial photonic 
magnetic fields due to their all-optical controllability. 
Moreover, their rich non-linear (quantum) dynamics [10] 
could be explored in the presence of an artificial mag¬ 
netic field. In general, the very flexible optical control 
could be used to create and explore novel features, e.g. 
varying the optomechanical coupling strength spatially 
and/or temporally, both adiabatically and with sudden 
quenches. Moreover, a second strong control laser could 
be used to create a spatially and temporarily varying op¬ 
tical on-site potential landscape. 
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Appendix A: Derivation of the effective magnetic 
Hamiltonian for the modulated link scheme 

Here, we derive the effective Hamiltonian 
hJ e ^e l( ^a} B dA + h.c. that describes the tunneling of 
photons from site A to site B in the presence of an 
effective magnetic field created using the modulated 
link scheme. We start from the full time dependent 
Hamiltonian Eq. <§• Since this second-quantized 
Hamiltonian is particle conserving we can switch to 
a first-quantized picture in the standard way. The 
corresponding single-particle Hamiltonian Hm reads 

(u A 0 ~J \ 

H m = ft 0 ljb -J 

\—J —J djj + 2<7q|/?| cos(Dt + (j))J 

It acts on the photon wavefunction \^) = (') 
where ip s describes the probability amplitude that the 
photon is localized on site s, s = A, B , /. Since the 
Hamiltonian is time periodic, there is a complete set of 
quasi-periodic solutions of the Schri'^oedinger equation, 
\^j(t + 27t/ 0)) = exp[— i2'KUjj/Yi]\^j{t)) where 2tt/YI is 
the period and index j spans the Hilbert space, j = 
1,2,3. In practice, one solves the eigenvalue problem 
£jm\4>jm)) = Hjm\^jm)) where H m -ihdt + Hm is 
the Floquet-Hamiltonian, £j m = h(ujj + mYt) are the 
quasienergies and \4>jm)) = exp [i(uj + mCl)t\\^j(t)) are 
time-periodic states, the so-called Floquet eigenstates 
[m £ Z] [45] . Notice that the Floquet Hamiltonian can be 
regarded as an operator on the extended Hilbert space of 
the time periodic vectors equipped with the scalar prod¬ 
uct 

((4*3 10m)) = t f (Al) 

In this framework, we can use the standard quan¬ 
tum mechanical perturbation theory to derive an effec- 
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UJj -\- U 



UJl 


Figure 5: Floquet level scheme, i.e. quasienergy levels, de¬ 
scribing the hopping transitions between sites A and B me¬ 
diated by the virtual transition trough site I accompanied by 
the exchange of a phonon. 


tive time independent single-particle Hamiltonian. We 
assume a resonant drive uoa ~ and weak tun- 

neling/driving, Ja,b/\^a,b ~ w/|, g 0 /3/£l <C 1. We iden- 
tify resonant Floquet-levels with quasienergies Huja and 
+ O) coupled via the third order virtual tunneling 
process through the interface site I shown in Figure [5| 
Up to leading order in perturbation theory, we can focus 
on the block of the Floquet Hamiltonian comprising the 
four unperturbed quasienergy levels that are involved in 
this process, cf. Figure [5j 


(UA 

0 

-J 

0 

0 

ujb A U 

0 

-J 

-J 

0 

UJI 

9o0 

V o 

-J 

9o0* 

OJl + 


Application of a standard Schrieffer-Wolff transformation 
[46ti48] . i.e. applying degenerate perturbation theory to 
third order, leads to the effective block diagonal Floquet 
Hamiltonian 


^ eff h (j eS A e i * w^+Wex) ‘ (A2) 

where Co s = uo s + J^/(uo s — ujj) with s = A, B and 
Jett = 9o\P\J aJb/[{u a -ui){u)b-ui)]. Finally, we turn 
back Hamiltonian | A2| into its second-quantized form and 
switch to a frame rotating with frequency Co a (&b) on 
site A (B). For a resonant drive, Co a = Cob + this 
yields the desired form of the second-quantized effective 
Hamiltonian, J e ff {e~ l ^a\aB + e^d^d^). 


Appendix B: Transmission amplitudes and density 
of states for the modulated link scheme 

Here, we calculate the LDOS for the modulated link 
scheme which is plotted in Figure [2] We use the full time 
dependent Hamiltonian Eq. <© extended to the whole 
lattice (including also the sublattice formed by the link 
sites). Since we are dealing with a time periodic system 
where the energy is not a constant of motion, we have to 
appropriately generalize the definition of the LDOS. A 
natural generalization of the standard definition to time- 
periodic systems is the following, 


where G(cj, m; j, 1) is the Floquet Green’s function 

pT poo 

G=^~ dr dte i( ' u+mn ' ,t+imnT {[dj(t + r),Oj(r)]). 
T Jo Jo 

The Floquet Green’s function describes the (linear) re¬ 
sponse of the array to a probe laser. More precisely, 
the light amplitude on site j in the presence of a probe 
drive on site 1 with frequency oo and amplitude a^ in ^ 
[described by the additional Hamiltonian term Hj = 
+ h.c.)\ is 


(fij(t)) = ^e-^ +m ^iv^^ n) G(cj,m;j,l). 


This is essentially a generalization of the Kubo formula 
which applies to any time periodic Hamiltonian. Using 
the input-output relations, a^ out \t) = a- zn \t) — y^dj, 
we can also calculate the field outside the cavity, 

(. d\ out \t )) = ^e- <( “ + "* n)t to(w,m;j, 1 )a (<n) , 


where 


to(u>, m; j, 1) = <5ji<5 m0 - inG{w, to; j, 1) 


is the transmission amplitude of a photon from site 1 
to site j if it has been up-converted m -times (or down- 
converted | m |-times for m negative). 

For a time-periodic system with a particle conserving 
Hamiltonian, the Floquet Green’s function can be easily 
expressed in terms of the first-quantized Floquet Hamil¬ 
tonian % = —idt — H(t), 

G(u>, to; j, 1) = <(j, to| (lj - H + iK/2)- 1 |l,0)) 


Notice that the Floquet Hamiltonian and the Green’s 
function can be regarded as operators acting on the ex¬ 
tended Hilbert space of the time -periodic photon states 
with the scalar product Eq. (A1 ). As such they acts on 
the time periodic states |j, m)), where index j indicates 
the lattice site and m the Fourier component. Thus, the 
density of states can be readily computed by diagonaliz¬ 
ing the Floquet Green’s function. We find 


pH = E z 

k V 


Wfc) 2 + « 2 /4 




where hoo^ are the quasienergies and \(/)k}} are the cor¬ 
responding Floquet eigenstates obtained by numerically 
diagonalizing H. Taking into account that the Floquet 
eigenfunctions | </>&)) forms a complete orthonormal basis 
of the Hilbert spac e of the time-periodic states [with the 
scalar product Eq. ( |A1| )], it immediately follows that the 
density of states is appropriately normalized, 


/ oo 

duop(po) = 27r. 

-oo 


p(u, j) = —2ImG(w,0;j,j) 










Appendix C: Transmission amplitudes for the 
frequency-conversion scheme 

For the frequency-conversion scheme we start from the 
linearized Langevin equations for the full array including 
the mechanical links modes hi hhi, 

k = itr 1 2 3 4 5 6 [H,b k \-rbk/2 + Vfi>£ n \ 

Oj = Oj] — KOj/2 + -y/KOj 111 ). (Cl) 

The first line (second line) describes the sites hosting a 
mechanical (optical) mode. The Hamiltonian H is given 
by Eq. 0 extended to the full array and the noise forces 
have the usual commutation relations |lj. Notice that 
Eq. is written in a frame where the optical modes on 
sublattice A and B are rotating with frequency ujli and 
CJL 2 , respectively. A probe laser on site 1 with frequency u) 
and amplitude a in is described by the additional Hamil¬ 
tonian term Hi = ih^a^ n \a\e~ lApt — h.c.), where 
A p = uj — ujls (s = 1, 2 for 1 on sublattice A or B, respec¬ 
tively). The linear response of the light amplitude on site 
j to such probe laser is given by the Kubo formula 

(«j (*)) = iV^cS m) e~ lIlpt Gaa t(Ap,j,l) 

-V^ (in) e iApt Gaa(-A p ,j,l), (C2) 

with the Green’s functions 

nOO 

Gaat(w,j,l) = -i / dte^da^t),^^ 0)]), 
poo 

Gaa(w,j,l) = ~i / dte wt {[a- } {t), ai(0)]). 

Jo 


Notice that in Figure 3 and 4 of the main text we plot the 
resonant p art o f the response corresponding to the first 
line of Eq. (C2). If j and 1 lie on different sublattices, the 
frequency of the prob e sig nal is converted [to read off this 
frequency from Eq. ( |C2| ), one has to keep in mind that 
the frame of reference is rotating at different frequencies 
on the two optical sublattices]. Finally, we note that 

'(outA 


the light transmitted outside of the sample 
t( j, 1) ( d[ m M can be readily computed using the input 


a • 


- WKCli 


output relations [49] dj = 

(C2) we find the transmission amplitude 


From Eq. 


to{u, 1, j) = <5ij - ikG m t(w,l,j). (C3) 


Since the transmission amplitudes of a probe laser beam 
are generally proportional to the corresponding light am¬ 
plitudes inside the array (on all sites except for the one 
where the light is injected), the amplitude patterns shown 
in Figures 3 and 4 could be directly measured by a posi¬ 
tion resolved measurement of the light scattered by the 
array. 

In order to calculate the transmission in Figures 3 and 
4 we have calculated the Green’s function numerically. 
We note that for an array with N x N optical sites, there 
is a total of TV(27V—1) sites (including also the mechanical 
sites) and a total of 27V(27V— 1) degrees of freedom. Thus, 
computing numerically the Green’s function amounts to 
inverting a 27V(27V — 1) x 27V(27V — 1) matrix. In Figure 
3 and 4 we have chosen TV = 22. 
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